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We study the motion of a particle in a periodic potential with Ohmic dissipation. In D = 1 
dimension it is well known that there are two phases depending on the dissipation: a localized phase 
with zero temperature mobility ^ = Q and a fully coherent phase with unaffected by the periodic 
potential. For D > 1, we find that this is also the case for a Bravais lattice. However, for non 
symmorphic lattices, such as the honeycomb lattice and its D dimensional generalization, there is 
a new intermediate phase with a universal mobility /i*. We study this intermediate fixed point in 
perturbatively accessible regimes. In addition, we relate this model to the Toulouse limit of the D + l 
channel Kondo problem. This mapping allows us to compute /i* exactly using results known from 
conformal field theory. Experimental implications are discussed for resonant tunneling in strongly 
coupled Coulomb blockade structures and for multi channel Luttinger liquids. 
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The quantum mechanics of a particle in a periodic po- 
tential coupled to a dissipative environment is a funda- 
mental problem in condensed matter physics |l| . A sim- 
ple theory based on the Caldeira-Leggett model of Ohmic 
dissipation was proposed in the mid 1980's as a possible 
description of the motion of a heavy charged particle in 
a metal [^j. In a one dimensional periodic potential it 
was shown that there are two zero temperature (T = 0) 
phases. For weak friction, the particle diffuses freely as 
if the periodic potential were absent. When the friction 
exceeds a critical value, however, the particle is localized 
in one of the minima of the potential. 

Recently there has been renewed interest in this quan- 
tum Brownian motion (QBM) model in connection with 
quantum impurity problems and boundary conformal 
field theory Q . It is isomorphic to the problem of tunnel- 
ing through a barrier in a Luttinger liquid, which is rele- 
vant to experiments in quantum wires |^ and tunneling 
in quantum Hall edge states |^ . Here the "coordinate" of 
the "particle" is the number of electrons that tunnel past 
the barrier. The periodic potential arises from the dis- 
creteness of the electron's charge. The Luttinger liquid's 
modes play the role of the dissipative bath. The particle's 
mobility corresponds to the electrical conductance. 

There are often multiple electron channels, due to spin 
and transverse degrees of freedom. The impurity problem 
then maps to a multi-dimensional periodic potential. In 
addition to the extended and localized phases, in two di- 
mensions it has been shown that there can be additional 
non trivial phases which may be accessed by tun- 

ing to a resonance. Using a similar analysis, Furusaki and 
Matveev recently found a similar intermediate phase in 
a model of resonant tunneling through a Coulomb block- 
ade structure They argued that the resonance fixed 
point is that of the multi channel Kondo problem. 

In this paper we consider the general problem of QBM 
on periodic lattices. We show that the lattice symme- 
try plays a crucial role in determining the T = phases. 
For the honeycomb lattice and its D dimensional gener- 



alization, there is a T = phase described by an inter- 
mediate fixed point, which we relate to the D -\- 1 chan- 
nel Kondo fixed point. Exploiting the mapping onto the 
Kondo problem, we compute exactly the fixed point mo- 
bility and critical exponents by borrowing results from 
conformal field theory. 

Integrating out the bath degrees of freedom, the QBM 
model is described by the Euclidean action 111 
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where r is the coordinate of the particle and vq, are di- 
mensionless Fourier components of the periodic potential, 
defined at the reciprocal lattice vectors G. (G is defined 
so that G • R is an integer for any lattice vector R.) The 
coupling to the dissipative bath is described by 
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where Tc is a short time cutoff. The friction is propor- 
tional to the coefficient of this term. However, by rescal- 
ing r and G this coefficient may be fixed. The lattice 
constant thus controls the strength of the friction. In 
accordance with Ref . 1^ , we define the dimensionless pa- 
rameter g = |R,„in|~^, where |Rmin| is the Bravais lattice 
constant, g is inversely proportional to the friction. A 
1 + 1 dimensional version of this theory has recently been 
analyzed by Kondev and Henley |l^ . 

Our system may be characterized by the mobility, 
which describes the average velocity of the particle in 
response to a uniform applied force. We define the di- 
mensionless mobility /i as the ratio of the mobility to the 
"perfect" mobility obtained in the absence of the periodic 
potential. When = 0, /i = 1. /i may be computed 
from linear response theory. 
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The effect of the periodic potential may be analyzed 
perturbatively in either of two limits. A weak potential 
may be studied by considering the renormalization group 
(RG) flows to leading order in vg, 



(a) 



dvG/di=il-\G\^)vG. 
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Clearly, if the shortest reciprocal lattice vector satisfies 
|Gniin| > Ij then all vg are irrelevant. The "small bar- 
rier" limit in which the particle diffuses freely is thus 
perturbatively stable. On the other hand, if |Gmin| < 1, 
then the system flows to a different strong coupling phase. 

When the barriers are large, the particle is localized in 
one of the minima of the potential with a small probabil- 
ity for tunneling to another. It is then more natural to 
consider a dual representation in which the the partition 
function is expanded in powers of the "fugacity" of these 
tunneling events For a Bravais lattice, this may be 

generated by expanding the dual action, 



^= R 
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tji may be interpreted as the matrix element for the par- 
ticle to tunnel between minima connected by a lattice 
vector R. Equivalently, k(T) describes the particle's tra- 
jectory in momentum space in a potential with the sym- 
metry of the reciprocal lattice. The RG flows to leading 
order in are then 



dtn/di = (1 - IRHtR. 



(6) 



The "large barrier" phase is thus perturbatively stable 
provided the shortest lattice vector satisfles |Rmin| > 1- 

For a one dimensional lattice, |RminP — I/5 and 
|GminP = 9- Thus either the small or the large bar- 
rier limit is stable, but not both. There are two phases: 
for g < 1 the system is localized and for g > 1 the sys- 
tem has perfect mobility. Clearly, this is also the case in 
higher dimensions for a lattice with cubic symmetry. 

In contrast, for a triangular lattice, IRminp = 
but IGminP = 4g/3. It follows that for 3/4 < 5 < 1, both 
the small and large barrier limits are stable. There must 
therefore be an unstable flxcd point separating the two 
phases, as indicated in Fig. la. A similar intermediate 
flxed point occurs in the single barrier problem of a spin 
1/2 Luttinger liquid (t). 

A perturbative analysis of this fixed point is possible 
for small v in the vicinity of 5 = 3/4. Specifically, con- 
sider a model with vg = v for the six "nearest neighbor" 
reciprocal lattice vectors. For v > 0, this produces a 
potential with minima forming a triangular lattice. For 
g — 3/4(1 -I- e), the RG flow to second order in v is 



dv/de =^ -ev + 2v^ 



(7) 



Provided v > and e > 0, there is an unstable flxed point 
V* = e/2, with RG eigenvalue e. The dimensionless mo- 
bility at this fixed point is universal, /i* = 1 — (37r^/2)e^. 

A similar analysis is possible in the dual theory for 
small t. Since the dual lattice is also triangular, the 
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FIG. 1. Flow diagrams for the (a) triangular lattice and 
(b) honeycomb lattice. The top (bottom) hnes represent the 
small V (t) limits. Stable (unstable) fixed points are de- 
picted by solid (dotted) lines, and arrows indicate the RG 
flows. The fixed point mobility is known perturbatively near 
g — 1/3, 3/4, 1 and at the points indicated by the dashed lines. 



results are identical, given the substitutions r ^ k, 
G — > R, g — !■ 3/(45) and v t. /i is mapped to 1 — fi. 
For g = 1 — e, the flxed point is at t* = e/2, with an 
exponent e, and mobility ^* = (37r^/2)e^. For g = v^/2, 
the theory is self dual, which implies that the fixed point 
mobility is /i* = 1/2. Piecing these results together, we 
obtain the fiow diagram in Fig. la. 

When V < 0, the minima of the potential described 
above form a honeycomb lattice. The honeycomb lat- 
tice is equivalent to the triangular lattice described above 
with a two site basis. If the triangular lattice constant 
is l/y^, then we again have |GminP = 4g/3. However 
|R„,i„|2 ^ l/(3.g) is now shorter. Thus, for 1/3 < g < 3/4 
both the large and small barrier limits are unstable, so 
that there must be a stable flxed point describing a new 
intermediate coupling phase. With v negative, Eq. (1) 
may also be viewed as the tight binding representation 
of a triangular lattice with tt flux per plaquette, whose 
reciprocal is the honeycomb lattice. Our analysis ap- 
plies to this dual theory as well. Callan et. al. |^ have 
recently found similar intermediate phases on a square 
lattice with magnetic flux. 

A perturbative analysis is again possible in the large 
and small barrier limits. For small barriers the fixed 
point of Eq. (7) is stable for u < and e < 0. The RG 
eigenvalue and mobility at the flxed point are the same 
as above. In the large barrier theory we must keep track 
of the two site basis of the honeycomb lattice. There are 
three nearest neighbors R for each site on the A sub- 
lattice. For the B sublattice the nearest neighbors are 
— R. The tunneling must alternate between the sublat- 
tices. This can be incorporated in the dual theory by 
introducing a spin 1/2 degree of freedom. For nearest 
neighbor hopping the dual action is then 



^= R 
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where R are among the 3 nearest neighbor lattice vectors 
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FIG. 2. Lattice of spin states for tlie 3 cfiannel Kondo prob- 
lem, wliicfi form two planes with constant Sf + 5*1 + 51 . 

of sublattice A and are spin 1/2 operators, Tlie 
intermediate fixed point may now be accessed perturba- 
tively for g — (1/3)(1 + e). We have computed the RG 
flow equation to order t^, 

dt/M = et-3t^. (9) 

For e > there is a stable fixed point at t* = -^/e/S, with 
RG eigenvalue 2e. The fixed point mobility is fj,* = tt^c. 

The flow diagram for the honeycomb lattice as a func- 
tion of g is summarized in Fig. lb. Unlike the cubic and 
triangular Bravais lattices, the T — mobility does not 
exhibit a discontinuous jump from to 1 as 5 is increased. 
Rather, the mobility interpolates smoothly between the 
two limits in the intermediate phase for 1/3 < g < 3/4. 
Below we show that for g — 1/2, the intermediate fixed 
point is that of the 3 channel Kondo problem. Exploit- 
ing this mapping, the mobility for g = 1/2 may be com- 
puted exactly. The generalization of the above analysis 
to other lattices and higher dimensions is straightforward 
and will be given elsewhere. In general, the existence of a 
stable intermediate phase requires a non symmorphic lat- 
tice symmetry, with a vector connecting equivalent sites 
that is shorter than any lattice translation. 

We now relate the stable intermediate fixed point to 
the multi channel Kondo problem by identifying the lat- 
tice symmetry in the Kondo problem. The Hamiltonian 
of the anisotropic TV channel Kondo model is , 

a.s i,a 

where a, s, i are channel, spin, and space indices, Sl^^ 
is the impurity spin, and sl{0) = i/'L(0)(f^L'/2)'/'as'(0) 
is the electronic spin in channel a at a; = 0. We con- 
sider an anisotropic model, characterized by dimension- 
less couplings Jz and = Jy = J±. Our analysis closely 
parallels that of Emery and Kivelson for the two channel 
Kondo problem |l^. We first bosonize the theory, and 
then do a rotation in spin space which transforms the Jz 
interaction. Upon integrating out the degrees of freedom 
away from x = 0, we obtain a theory in terms of the 
boson fields at the impurity which closely resembles the 
lattice models studied in this paper. The details of this 
mapping will be presented in a longer article, however its 
essence may be understood quite simply. 

When J± = the states of the system may be charac- 
terized by the total spin in each of the N channels. 



The possible values of form a A'' dimensional cubic 
lattice. J_L "hops" the system between sites on this lat- 
tice. Since H conserves the total spin of the electrons 
plus the impurity, the system is constrained to lie on one 
of two lattice planes with constant S^^^ + where 
S'i^p = ±1/2. For N = 3 each lattice plane forms a 
triangular lattice, as sketched in Fig. 2. Viewed from 
the (111) direction, the two lattice planes form a "cor- 
rugated" honeycomb lattice in which the two triangular 
sublattices are displaced in the perpendicular direction. 
For general TV, the "lattice planes" consist of two inter- 
penetrating iV — 1 dimensional close packed lattices. For 
N = 4, they form a corrugated diamond lattice. 
Now consider QBM on such a lattice described by 

^ = S'o[k]- [ Jr+e^2.(R|,k||+j^,fe,)^/^^^; 
J ''c „ L 

R|| 

For N = 3, k is a 3 dimensional vector with compo- 
nents k|| and k± parallel and perpendicular to the lattice 
plane. R|| are chosen from the 3 nearest neighbor lat- 
tice vectors for the honeycomb lattice, and R± is the 
perpendicular displacement between the two sublattices. 

guarantee that the hopping alternates between the 
two lattice planes. This model is identical to the multi 
channel Kondo problem, with t — J_l/2, provided the 
lattice constants are chosen to give the appropriate scal- 
ing for J±. For Jz = 0, the dimension of s+(0) is 1, so 
the cubic lattice constant in Fig. 2 is unity. It follows 
that |R|| p = 1 — 1/N. If the lattice constant of the close 
packed Bravais lattice is l/y/g, then g = 1/2. Finite 
Jz may be treated nonperturbatively using bosonization 
]l2[ , and affects the dimension of s+(0). This leads to 
a distortion the lattice in the perpendicular direction, 
Rl = (l- NJzl2f/N. Note that = for Jz = 2/N, 
so that the perpendicular direction decouples. This is the 
N channel generalization of the Toulouse limit [p^-p^. 
A central point of this paper is that this limit of the N 
channel Kondo model is identical to the g = 1/2 QBM 
model on a A^ — 1 dimensional "honeycomb" lattice. 

Note that the motion perpendicular to the planes alter- 
nates, whereas the motion parallel to the planes does not. 
This gives rise to a renormalization of (or equivalently 
Jz) but not R||. This may be seen from a RG analysis 
similar to that of Anderson, Yuval, and Hamann ||l5|,|l6| . 
Expressing the flow equations in terms of Jz and , we 
find to order J^, 

dJz/di = Ji [1 - (A^/2) J,] , (10) 
dJ^/di = J^J, [1 - (A/4) J,] - (Af/4) Ji. (11) 

The RG flows are shown in Fig. 3. Jz flows towards 
2/A^, the Toulouse limit, shown by the dashed line. For 
A = 3 this is the same as the dashed line in Fig. lb. The 
intermediate fixed point for the honeycomb lattice with 
g — 1/2 is the same as the fixed point of the 3 channel 
Kondo problem. Varying g adiabatically connects the 
multi channel Kondo fixed point to the strong and weak 
barrier limits described perturbatively in (7) and (9). For 
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FIG. 3. Flow diagram for the A'^ channel Kondo model for 
small Jx. The dashed line is the Toulouse limit, Jz = ijN . 
The strong coupling fixed point is marked with the full circle. 



large TV, the fixed point dX — = 2/N approaches 
the strong barrier limit, and is perturbatively accessible 
p^ . Perturbations which break the symmetry between 
the two sublattices are relevant and act like a magnetic 
field in the Kondo problem. The system then flows to a 
fixed point of the lower symmetry lattice. 

Conformal field theory allows for an exact description 
of the multi channel Kondo fixed point We can 

thus identify the critical exponents and mobility for the 

— 1 dimensional generalized honeycomb lattice model 
for g = 1/2. The RG eigenvalue of the leading irrel- 
evant operator at the fixed point is —2/{N + 2). The 
mobility is computed by identifying the appropriate cor- 
relation function in the Kondo model. The analogue 
of r in the Kondo model is the spin in each channel, 
S^, projected onto the TV — 1 dimensional lattice plane 
with constant. This corresponds to the operator 

O = '^ls'^ss''^aa''4^a's'> where T"^ is one of the TV - 1 diag- 
onal generators of SU{N). The "current" r corresponds 
to a flow of spin between the different channels. Ludwig 
and Affleck have computed all correlation functions 
of O exactly. Borrowing their results, we obtain for the 
fixed point mobility. 



/X* = 2 sin^ 



TV- 



(12) 



For TV = 2, the Kondo fixed point has fi* — 1, and is at 
the "small barrier" limit. For TV = 3, this value is plotted 
in Fig. lb. For large TV, the mobility and RG eigenvalue 
may be found perturbatively in a manner analogous to 
the e expansion following Eq. (9). We have checked 
that they agree with the exact result to leading and sub 
leading order in 1/TV. 

Recently, Furusaki and Matveev |^ have studied 
Coulomb blockade resonances in a spin degenerate quan- 
tum dot with quantum point contact leads. For A <C 
T ^ Ec, (A is the dot's level spacing and Ec is the 
Coulomb charging energy), they argue that resonances 
are controlled by the 4 channel Kondo fixed point. This 
mapping may be understood in terms of the lattice of 
allowed charge states for the dot and the four lead/spin 
channels: a symmetric dot on resonance has the sym- 
metry of a diamond lattice with g = 1/2. Our analysis 
allows us to identify the universal on resonance conduc- 
tance. A voltage between the leads corresponds to a force 
F = eV, and the resulting current is / = e|r|. From (12), 
with TV = 4, /i* = 1/2, leading to G* = {l/2)e^/h M. 



Our analysis also applies to resonant tunneling through 
a single resonant state (i.e. T < A), for a multi channel 
Luttinger liquid with repulsive interactions 0|. For two 
channels, a tunneling barrier maps to a two dimensional 
lattice with g < 1, which, when tuned to a resonance, 
can in principle have a non symmorphic distorted hon- 
eycomb lattice symmetry. Such resonances would be an 
adiabatic cousin of the 3 channel Kondo effect. Analo- 
gous resonances have been observed in a single channel 
Luttinger liquid ||^. 

In summary, we have presented a general theory of 
quantum Brownian motion on periodic lattices. For the 
honeycomb lattice and its D dimensional generalization, 
there is a non trivial intermediate phase, which we have 
identified with the multi channel Kondo problem. Pre- 
sumably, other non symmorphic lattices also display such 
phases, and it would be interesting to classify them using 
conformal field theory. 

We thank M.P.A. Fisher, K.A. Matveev and E.J. Mele 
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supported by NSF grant DMR 95-05425. 
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